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Abstract —In this article, we propose an improved mnitiple 
feedback successive interference cancellation (IMF-SIC) algo¬ 
rithm for symbol vector detection in multiple-input multiple- 
output (MIMO) spatial multiplexing systems. The multiple feed¬ 
back (MF) strategy in successive interference cancellation (SIC) 
is based on the concept of shadow area constraint (SAC) where, 
if the decision falls in the shadow region multiple neighboring 
constellation points will be nsed in the decision feedback loop fol¬ 
lowed by the conventional SIC. The best candidate symbol from 
multiple neighboring symbols is selected using the maximum 
likelihood (ML) criteria. However, while deciding the best symbol 
from multiple neighboring symbols, the SAC condition may occur 
in subsequent layers which results in inaccurate decision. In 
order to overcome this limitation, in the proposed algorithm, 
SAC criteria is checked recursively for each layer. This results 
in successful mitigation of error propagation thus significantly 
improving the bit error rate (BER) performance. Fnrther, we 
also propose an ordered IMF-SIC (OIMF-SIC) where we use log 
likelihood ratio (LLR) based dynamic ordering of the detection 
sequence. In OIMF-SIC, we use the term dynamic ordering in the 
sense that the detection order is updated after every successful 
decision. Simulation resnlts show that the proposed algorithms 
outperform the existing detectors such as conventional SIC and 
MF-SIC in terms of BER, and achieves a near ML performance. 

Index Terms —MIMO detection, spatial multiplexing, decision 
feedback, successive interference cancellation, bit error rate. 

I. Introduction 

With the rapid growth in demand for higher data rates, wire¬ 
less technology is shifting towards the systems with multiple 
antennas. Multiple-input multiple-output (MIMO) systems in 
wireless communications provide significant improvements in 
wireless link reliability and the achievable capacity Q-S- 
Recently, the systems with large number of antennas (tens 
to hundreds of antennas) called as large-MIMO systems are 
getting increased attention Q because they provide additional 
multiplexing and diversity gains. Using spatial multiplexing 
in MIMO systems, multiple data streams can be transmitted 
simultaneously from different transmit antennas. However, the 
detection of these data streams at the receiver end is challeng¬ 
ing. To achieve minimum bit error rate (BER) performance 
in MIMO systems, an exhaustive search over all the possible 
transmit vectors is required. This is called as the maximum 
likelihood (ML) search for symbol vector detection in MIMO 
systems. However, as the number of antennas grow, the num¬ 


ber of possible transmit vectors also increases exponentially 
and thus the ML search becomes computationally impractical. 
Sphere decoder, a well known MIMO detector achieves near 
ML performance but is practical only up to limited number of 
dimensions 0. 

The traditional MIMO detection techniques involve lin¬ 
ear detection such as zero forcing (ZE) detector and mini¬ 
mum mean squared error (MMSE) detector. These detectors 
use linear transformation of the received vector in order to 
estimate the transmitted symbol vector. Comparatively the 
MMSE detector is superior over ZE detector in terms of 
BER performance, but still their performance is far inferior 
compared to the ML performance. Another type of detection 
involve non-linear detectors such as vertical Bell labs layered 
space time architecture (V-BLAST) ||^ which utilizes ordered 
SIC for symbol vector detection. Some of the other MIMO 
detection technique involve ant colony optimization based 
MIMO detection Q, message passing based algorithm Q, 
channel hardening based algorithm Q, and lattice reduction 
aided MIMO detection i fTOl , | fTT| . As an alternative, in this 
study we employ SIC based detection technique for MIMO 
detection problem. 

SIC based MIMO detection is a well known detection 
technique where symbols are detected sequentially |12|. In 
SIC, after detecting each symbol, its interference is canceled 
from the received vector in order to improve the instantaneous 
signal to interference plus noise ratio (SINR) for the remaining 
symbols. It is also known as layered detection where in 
each layer, one symbol is detected. However, it suffers from 
error propagation which occur due to wrong decisions in 
early stages of the algorithm p3| , p^ . Eor improving the 
performance of SIC, a well known technique known as V- 
BLAST 1^ technique is used which performs the following 
steps: (1) SNR based ordering in the detection sequence, (2) 
using ZE or MMSE for nulling the interference among the 
data streams and (3) detecting a symbol and canceling its 
interference i.e. SIC. Eurther, to mitigate the effect of error 
propagation, one of the available algorithms in the literature 
include multi-branch SIC (MB-SIC) algorithm |15| where 
the concept of multiple branch (MB) processing is utilized. 
Each branch of MB-SIC differs in ordering of the detection 
sequence thus resulting in a higher detection diversity over the 




conventional SIC. Recently, a multiple feedback (MF) strategy 
is proposed in | fT^ for SIC based MIMO detection. Multiple 
feedback SIC (MF-SIC) algorithm proposed in is based 
on the concept of shadow are constraint (SAC) where if the 
decision fall in the shadow region (unreliability region) then 
multiple constellation points are used in the decision feedback 
loop. The best symbol from the multiple symbols used in the 
decision feedback is selected using the ML criteria. 

In this article, we propose an improved multiple feedback 
successive interference cancellation (IMF-SIC) for symbol 
vector detection in MIMO spatial multiplexing systems. In 
MF-SIC, once a decision falls into the shadow region, conven¬ 
tional SIC is used to find the best one symbol for the corre¬ 
sponding layer. However, it may happen that while computing 
the best symbol for a given layer, multiple decisions from the 
subsequent layers fall into the shadow region. This condition 
is not checked in MF-SIC which sometimes results in error 
propagation and limits the BER performance. In improved 
MF-SIC (IMF-SIC), we overcome this limitation by checking 
the shadow region criteria recursively rather than using the 
conventional SIC for searching the best candidate symbol. This 
results in significant reduction in error propagation and thus 
the BER performance can be improved. Also, we propose an 
ordered IME-SIC (OIME-SIC) algorithm where we employ 
the log likelihood ratio (LLR) based dynamic ordering in the 
detection sequence GZl- By the term dynamic ordering we 
mean that after every successful decision about a symbol, 
the ordering of the detection sequence is updated based on 
new LLR values. In OIME-SIC, same ordering pattern is 
followed even for searching the best symbol when a decision 
is found unreliable. Simulation results show that the proposed 
algorithms significantly outperform the conventional SIC and 
the ME-SIC based detection methods, and that they achieve a 
near ML performance. 

Rest of the article is organized as follows. In Sect. |II] we 
present mathematical model of MIMO system. An overview 
of some of the traditional MIMO detection methods is given in 
Sect. Ill The proposed algorithms for symbol vector detection 


in MIMO spatial multiplexing systems is discussed in Sect. IV 


In Sect. |V] we show the simulation results on comparison of 
BER performance. Einally, in Sect. VI we conclude the article. 



Fig. 1. MIMO system model. 


hi j for i = 1,2,..., Nj. and j = 1,2,..., Nt assumed to 
be independent and identically distributed (i.i.d.) as complex 
normal with mean 0 and variance 1 i.e. ~ CA/^(0,1). The 
received symbol vector y, after demodulation and matched 
filtering can be written as 

y = Hs-fn, (1) 

where n is the additive white Gaussian noise (AWGN) with 
its element Ui for i = 1,2,..., as i.i.d. ^ CN(Q,<j‘^). 
The signal to noise ratio (SNR) is defined as lOlogj^p ^ 
where Eg is the average energy per symbol defined as 
NtEg = tr(E[ss^]) and cr^ is the element-wise noise variance. 
tr(-) denote the trace of a matrix, E[-] denote the expectation 
operation and (•)^ denote the Hermitian of a matrix. The 
channel state information (CSI) is assumed to be perfectly 
known at the receiver but is unknown at the transmitter. The 
ML solution for a given received vector y and a known channel 
matrix H can be found as 

Sail = arg min ||y - Hs||^. (2) 

sGA”t 

ML detection (MLD) requires an exhaustive search over 
possible transmit symbol vectors where M is cardinality of 
the constellation set A. This search grows exponentially with 
increase in the modulation order (M) or increase in the number 
of transmit antennas (Nt) or both. 


IT System Model 

In this section, we discuss the mathematical model of 
MIMO system as shown in Eig[T] We consider a point to 
point spatially multiplexed MIMO link with Nt transmit 
antennas and N^. receive antennas. In spatial multiplexing, 
each antenna transmits a different information symbol Si for 
i = 1,2,... ,Nt . Let s = [si, S 2 , ■ • • , be the Nt x 1 
dimensional transmit vector, denote the transpose of a 
matrix. Each entry Si of the transmit vector s is taken from 
a signal constellation A, for example A = {—1 — 1*, —l-f 
li, 1 — li, 1 -f li} for 4-QAM signaling. The channel over 
which the symbol vector s is transmitted is assumed to be 
a Rayleigh distributed flat fading channel. Let us consider 
an Nr X Nt complex channel matrix H with its elements 


III. Overview of traditional MIMO detection 

TECHNIQUES 

In this section, we discuss some of the traditional MIMO 
detection techniques such as linear detection techniques which 
involve zero forcing (ZE) based MIMO detection and min¬ 
imum mean squared error (MMSE) based MIMO detection, 
and nonlinear detection techniques such as SIC based MIMO 
detection. 

A. Zero Forcing based MIMO detection 

In ZE based MIMO detection, a linear transformation matrix 
WzF is used in such a way that it removes the effect 
due to the channel impairments on the received vector y 
i.e. WzF = which is the pseudo inverse of 











channel matrix H. The linear transformation of the received 
vector using 'Wzf can be written as 

y = W^i^y = s + Wzfh, (3) 

Thus the ZF estimate of the transmitted symbol can be found 
as 

SZF = Q[W^Fy], (4) 

where Q[-] is the quantization operation which maps the soft 
values on to the nearest constellation point. The computational 
complexity of ZF receiver is very less, however the BER 
performance is far inferior to the ML performance. The main 
drawback of ZF based MIMO detection is the problem of noise 
enhancement which degrades its performance. 

B. MMSE based MIMO detection 

MMSE based MIMO detection is based on hnding a 
suitable linear transformation matrix which minimizes the 
mean squared error between the transmit vector x and the 
transformed received vector as 

Wmmsb = minE[||x - Wyf ]. (5) 

w 

The solution to Eq. can be written as 

^MMSE = (H^H + (6) 

Using the transformation matrix "Wmmse, the MMSE solu¬ 
tion can then be computed as 

SmMSE = Q[WMMSEy]- (7) 

MMSE based MIMO detection overcomes the problem of 
noise enhancement in ZF detectors. Thus it achieves better 
BER performance over ZF but still the performance of MMSE 
receiver is inferior as compared to the ML performance. One 
more drawback is that the MMSE solution needs knowledge 
of the noise variance cr^. 

C. SIC based MIMO detection 

In SIC based MIMO detection, the symbols transmitted 
from different transmit antennas are detected in a sequential 
manner. After detecting a symbol for a particular transmit 
antenna, its interference is canceled from the received vector in 
order to reduce the interference between different data streams 
at the receiver. The received vector in Eq. [T] can be rewritten 
as 

y = (hisi + h 2 S 2 H-h hAr^SArJ-f n. (8) 

Using the knowledge of channel matrix H a hlter is dehned 

2 

w, = (9) 

where H,_i is the matrix with its column taken from the 
i to Nt columns of the channel matrix H. We use vector 
Yi as the vector left after removing the interference due to 

■ ■ ■ ; ^2 — 1 

i-1 

y2 = y-^h,s,- (10) 

i=i 


Using Wi, the hrst symbol Si is estimated as 

Sz = Q[2)] = Q[wf y.,], (11) 

where Si is estimate of the transmitted symbol s^. w^, 
and Yi are updated after every successful decision about a 
symbol. SIC is an iterative detection technique, where in every 
iteration one symbol is detected. Thus a total of Nt iterations 
are performed to detect the complete symbol vector. A pseudo 
code of SIC based MIMO detection is given in Algorithm 

Algorithm 1 SIC based MIMO detection 
1 : input: y, H, Nt,Nr; 

2 : initialize: yo = y, Hq = H ; 

3: for z = 1 : 1 : iVt do 
4: Compute Wi = 

5: 2) = wfYi-l 

6: Sj = Q\zi] 

1: Update yi = yi_i - hiSi 

8: Update Hi = [hi+ihj +2 • • • HatJ 

9: end for 

10: output: s = [.S 1 S 2 • • • SatJ is output solution vector 


However, SIC suffers from the problem of error propagation 
which occurs due to the wrong decisions in early stage of 
the algorithm. Once an error occurs, it propagates to the later 
stages and thus increases the number of errors. This degrades 
the BER performance of SIC based MIMO detection. 

IV. Proposed algorithm for MIMO detection 

In this section, hrst we discuss the MF-SIC based algorithm 
for symbol vector detection in MIMO systems. We will then 
discuss the IMF-SIC algorithm and OIMF-SIC algorithm. 

A. MF-SIC algorithm 

The conventional SIC suffers from the problem of error 
propagation which degrades the BER performance. As a so¬ 
lution to mitigate the effect of error propagation, MF strategy 
has been proposed in p6) for SIC based MIMO detection. 
The MF strategy is based on the concept of reliability of a 
soft decision in SIC. Let us consider the soft decision 
within SIC as 

Zk = wf yfe. (12) 

The constellation point near to the soft decision Zk is nothing 
but its quantized value i.e. Q[2fe]. Let us dehne the distance 
between the soft decision and its quantized value as 

d=\zk-Q[zk]\ (13) 

If the soft decision in SIC is within a predehned threshold 
radius dth (fe. if d < dth) around the nearest constellation 
point then the decision is said to be reliable and its quantized 
value is used in the decision feedback. But if the decision is 
not within the radius (if d > dth) then multiple neighboring 
(say S) constellation points are used instead of the quantized 
value followed by the conventional SIC. It can be thought of 
as running parallel streams one for each symbol and at last 






the best one is selected using the ML metric (Eq. |^. MF-SIC 
performs superior than the conventional SIC in terms of BER. 
A pseudo code of MF-SIC is shown in Algorithm 


Algorithm 2 MF-SIC based MIMO detection 
1 : input: y, H, N-t, Nj., dth, S', 

2 : initialize: yo = y, Hg = H ; 

3: for j = 1 : 1 : do 

4: Compute Wj = -f 

5: = wf yj_i 

6: d=\Zi-Q[Zi\\ 

7: if (i < dth then 

8: Si = Q[Zi] 

9: else 

10 : for j = 1 : 1 : S' do 

11: = N^^\zi) i.e. initialize by neighbor¬ 

hood of Zi 
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for fc = i -f 1 : 1 : do 

Zk = wf yfc_i 
Sk = Q\zk\ 

Update Yk = Yk-i - ^kSk 
Update Hfc = \\ik+i^k +2 ■ ■ ■ h^v*] 

end for 

Compute the solution vector 

end for 


jopt = argminj=i,..,,s ||y - p 

s. = 

end if 

Update Yi =Yi-i- h^Si 
Update Hi = [hj+ih ,+2 • • • h^vj 

end for 

output: s = [siS 2 • • ■ satJ is output solution vector 


B. Proposed IMF-SIC algorithm 

In IMF-SIC, we propose an improved MF strategy for SIC 
based symbol vector detection in MIMO spatial multiplexing 
systems. In MF-SIC algorithm p6) , multiple feedback loops 
are followed by the conventional SIC to check for the optimal 
decision for a particular layer. However, while searching 
for the optimal decision, the Shadow region (uru'eliability) 
criteria is not considered in the lower layers and hence 
there are chances of selecting a sub-optimal decision about 
a symbol which degrades the BER performance. In order 
to improve the accuracy about a decision among multiple 
symbols in the feedback loop, we check the shadow region 
criteria recursively in IMF-SIC algorithm. We therefore 
formulate a separate sub-routine (function) named as imf sic 
in Algorithm which checks the unreliability condition and 
optimally decide the symbol. Sub-routine imf sic is called 
from the main function which is given in Algorithm and 
also from the sub-routine itself. We denote the number of 
calls of the sub-routine by L in IMF-SIC algorithm which 
is nothing but the number of times uru'eliability condition is 
checked once a decision falls into the shadow region. The 
number of neighbor'ing symbol which are used in the decision 


feedback is denoted by S. 


Algorithm 3 IMF-SIC for MIMO detection 
1: input: y, H, Nt, Nr, dth, S', 

2 : initialize: yo = y, Hg = H ; 

3: for j = 1 : 1 : A^t do 
4: Compute Wi = -f 

5: Zi = wf yj_i 

6: d=\Zi-Q[zi\\ 

1: it d < dth then 

8: Si = Q[zi] 

9: else 

10: Si = imf_sic(Hi_i,dt/,,yi_i,j) 

11: end if 

12: Update yi = yi_i - hiSi 

13: Update Hi = [hi+ih.j +2 • • • hjvJ 

14: end for 

15: output: s = [siS 2 • • • SNt] is output solution vector 


Algorithm 4 imf_sic 

1: function [si]=imf_sic(Hi_i, yi_i, i) 

2 : for j = 1 : 1 : S' do 

3: xp’* = Af^^^Zi) i.e. initialize by neighborhood of 

5i 
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for A: = z + 1 : 1 : A^t do 

if d < dth then 

Si = Q\Zi\ 

else 

Si = imf_sic(Hfc_i,dt/i,yfc_i,fc) 

end if 

Update yfc = Yk-i - hfeSfc 
Update Hfc = [hfc+ihfc+2 • • • h^vj 

end for 

Compute the solution vector x^A) 

end for 

jopt = argminj=i,...,s ||y - Hx(^)|p 

S. _ rr^iopt) 

Si — J-i 

return Si 


Significant reduction in the effect of error propagation is 
noticed using IMF-SIC (discussed in Sect. |V]) which results 
in improvement in the BER performance over SIC and the 
MF-SIC. However, IMF-SIC does not include ordering in the 
detection sequence which plays a very important role in mit¬ 
igating etTor propagation in SIC based detection techniques. 
In the next section, we discuss ordered IMF-SIC (OIMF-SIC) 
algorithm for MIMO detection. 

C. Proposed Ordered IMF-SIC (OIMF-SIC) algorithm 
In this section, we propose an ordered IMF-SIC (OIMF- 
SIC) algorithm which incorporates the LLR based ordering 
in the detection sequence. Error propagation can be mitigated 
successfully if the decisions in the early layers of SIC are 
detected correctly. In several SIC based detection algorithms. 
















channel norm ordering is used such as in V-BLAST technique 
. However, channel norm ordering does not consider the in¬ 
stantaneous noise in the received vector which has a signihcant 
impact on the BER performance. Thus a log likelihood ratio 
based ordering for multi-user detection in CDMA systems GZ) 
has been proposed in for MF-SIC based MIMO detection 
to overcome the sub-optimal behavior of channel norm based 
ordering. We employ LLR based dynamic ordering in the 
IMF-SIC algorithm to order the detection sequence. The term 
dynamic ordering is used in the sense that the ordering is 
updated after every successful symbol decision. However, the 
modihcation in the proposed method is that the ordering is 
also used in the oimf sic sub-routine which help in a more 
precise decision. The LLR value of a soft estimate Zj can be 
written as m 

L, = {l-hfR-%)-^\z,\ (14) 

where R = HH^ + a^I. The value of the argument j which 
maximizes Lj is arranged in descending order. Let us consider 
a set B be the set of indices of the symbols which are not 
detected so far. The set B is updated after every successful 
decision. For ordering the detection sequence in every iteration 
we consider the indices which are in the set B i.e. in Eq. 
j G B. Let us also consider the ordered sequence t which 
is updated in every iteration. In Algorithm]^ we present the 
OIMF-SIC algorithm for MIMO detection. 



Fig. 2. BER performance comparison for uncoded 4x4 V-BLAST MIMO 
system with 4-QAM and the value used for dti^ = 0.2, L = 2 and 3 = 4. 


Algorithm 5 OIMF-SIC for MIMO detection 


input: y, H, Nt, Nr-,dth, S; 

initialize: yi = y, Hi = H, R and B = {1,2,..., Nt} 
for TO = 1 : 1 : do 

for ( = 1:1: length{B) do 
Compute Wb, = (Hb,H^ + 

Zbi = ^bjh 


Compute Lb, using the equation 14 

end for 


Arrange the values Lb Nbj G B in descending order 
Update the detection order t 

Start detection using IMF-SIC given in Algorithm 
|3]vith ordered sequence t. 

12 : end for 

13: output: s = [siS 2 • • • Sat,] is output solution vector 


V. Simulation results 

In this section, we discuss the BER performance plots of 
the proposed algorithms and compare with the performance 
of the conventional SIC and the MF-SIC algorithm for 4x4, 
8x8 and 16 x 16 MIMO systems with 4-QAM and 16-QAM 
signaling. The results are generated using MATLAB and the 
BER is averaged over 10"^ Monte-Carlo simulations. 

In Fig. [2] we compare the BER performance of the con¬ 
ventional SIC and MF-SIC with the proposed algorithms for 
4x4 MIMO system with 4-QAM. The value of threshold 
parameter dth used is 0.2 and the number of recursions 
L = 2. Observation reveals that the proposed algorithms 


achieve a near ML performance and also the performance 
of IMF-SIC and OIMF-SIC algorithms is similar in terms 
of BER performance. At BER of 10“^, IMF-SIC results in 
approximately 1 dB gain in SNR over MF-SIC. In Fig. [3] the 
BER performance is compared for 8 x 8 MIMO system with 
4-QAM. A signihcant improvement in the BER performance 
of IMF-SIC and OIMF-SIC can be seen over SIC and MF- 
SIC based MIMO detection techniques. An SNR gain of 
around 2 dB can be achieved by using IMF-SIC over MF- 
SIC to achieve a target BER 2 x 10“^. Due to near optimal 
performance of IMF-SIC for 4 x 4 and 8x8 MIMO systems 
with 4-QAM the advantage of dynamic ordering cannot be 
observed. However, with increase in number of antennas or 
the modulation order, ordering plays an important role in 
mitigating the error propagation. 

Fig. [^presents the BER performance comparison for 16x 16 
V-BLAST MIMO system with 4-QAM. A signihcant improve¬ 
ment in the performance of OIMF-SIC can be seen when 
compared with IMF-SIC for L = 2. However, to achieve the 
near ML performance the number of recursions are increased 
to L = 3 and dtt = 0.5 for OIMF-SIC. Thus, the effect of 
increase in number of recursions on the BER performance can 
be observed. 

The BER performance for 4 x 4 and 8x8 V-BLAST 
MIMO systems with 16-QAM signaling is shown in Figs. 
and|^ In 16-QAM the total number of constellation points is 
16 and hence in order to get a reliable decision, we keep the 
number of neighboring contellation point in the feedback as 














Fig. 3. BER performance comparison for uncoded 8x8 V-BLAST MIMO 
system with 4-QAM and the value used for dth = 0.2, L = 2 and S = 4. 



Fig. 4. BER performance comparison for uncoded 16 X 16 V-BLAST MIMO 
system with 4-QAM. 



Fig. 5. BER performance comparison for uncoded 4x4 V-BLAST MIMO 
system with i6-QAM and the value used for = 0.2, L = 2 and S = 8. 


S = 8. The threshold value dth is kept same i.e. dth = 0.2. 
It can be observed that the BER performance of OIMF-SIC is 
close to within the 0.1 dB of the ML performance. Ordering 
in IMF-SIC results in almost 1 dB gain in SNR over the IMF- 
SIC without ordering. However, in Fig. to achieve the BER 
close to ML, the number of recursions is increased to L = 3 
in OIMF-SIC. This shows a SNR gain of more than 2 dB over 
IMF-SIC at a target BER 10“^. 

VI. Conclusions 

We proposed an improved multiple feedback successive in¬ 
terference cancellation (IMF-SIC) algorithm for symbol vector 
detection in MIMO systems. We also proposed an ordered 
IMF-SIC (OIMF-SIC) algorithm which employs log likelihood 
ratio (LLR) based dynamic ordering in the detection sequence. 
The ordering of the detection sequence is updated after every 
successful decision and thus the effect of error propagation 
is mitigated. The proposed algorithms significantly improves 
the BER performance and outperform the conventional SIC 
algorithm and the MF-SIC algorithm, and achieve near ML 
performance. 
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